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CN ■ Abstract 

In this paper we deal with diffusive relaxation limits of nonlinear systems of Euler 
type modeling chemotactic movement of cells toward Keller-Segel type systems. The 
^^ I approximating systems are either hyperbolic-parabolic or hyperbolic-elliptic. They all 

Cn ■ feature a nonlinear pressure term arising from a volume filling effect which takes into 

account the fact that cells do not interpenetrate. The main convergence result relies 
p, ' on compensated compactness tools and is obtained for large initial data under suitable 

.^ I assumptions on the approximating solutions. In order to justify such assumptions, we 

also prove an existence result for initial data which are small perturbation of a constant 
state. Such result is proven via classical Friedrichs's symmetrization and linearization. 
C^ ■ In order to simplify the coverage, we restrict to the two-dimensional case with periodical 

boundary conditions. 
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^ , 1 Introduction 

lO . This paper deals with diffusive relaxation limits for the nonlinear hyperbolic model describing 

chemotactic movement of cells, also known as the persistence and chemotaxis model. 



OO 



o 

oo 
O 



'drP + V-{pv) = 

drV + v-Vv + Vg{p) =Vc-dv (1.1) 

adrC = Ac + ap — j3c, 



^ . with a, /?, d, a positive constants. The function gi^p) is taken of the form gi^p^ = p'^ with 7 > 0, 

H \ we shall discuss this choice later on in this introduction. The model (jl.ip has been introduced 

and motivated very precisely in [?], whereas similar models have been also discussed in [?, ?, 
?, ?]. We shall briefly summarize the biological motivations behind (jl.ip by framing them in 
the general context of PDE systems describing chemotactical phenomena. 

The analysis of partial differential equations modeling chemotaxis goes back to the work of 
Keller and Segel [?], who proposed a macroscopic model for aggregation of cellular slime molds, 
and to the earlier related work of Patlak [?] , who derived similar models with applications to 
the study of long-chain polymers. In the successive decades, the term chemotaxis has been 
used to represent the dynamics of several biological systems (such as the bacteria Escherichia 
Coli, or the amoebae Dyctiostelium Discoideum, or endothelial cells of the human body 
responding to angiogenic factors secreted by a tumor) in which the motion of a species is 



biased by the gradient of a certain chemical substance. Typicahy these models consist of a 
system of drift -diffusion equations of the form 

(pt = Ap-V-ipxip,c)Vc) ^^2^ 

\ct = Ac + r{p,c), 

with diffusion terms modeling random motion for the density p of the individuals (cells, 
bacteria and so on) and for the concentration c of the chemoattractant (the chemical substance 
responsible of the chemotactical movement), first order drift terms modeling chemotactical 
aggregation and zero order reaction terms in the equation for the chemoattractant. The 
coefficient xip^ c) is called chemotactical sensitivity. The simplest model combining diffusion 
and chemotaxis is the well known parabolic-elliptic Patlak-Keller-Segel system (or simply 
Keller-Segel system) 

fpi = Ap-V-(xpVc) 

\o = Ac + p, 

which has the interesting mathematical feature of producing smooth solutions for small ini- 
tial norms (in the appropriate space) and blow-up (in the form of concentration to deltas) 
for large initial norms. The rigorous analysis of such mathematical issue (also extended to 
fully parabolic systems and to more complex models) has attracted the attention of many 
mathematicians in the last decades. We mention the pioneering works of Jager-Luckaus [?], 
Nagai [?], Herrero- Velazquez [?] among others. The existence vs. blow-up problem in two 
space dimensions for the classical Keller-Segel model (|1.3p has been completely solved in [?], 
where the authors proved that if the initial mass is less than a threshold value m* (depending 
on the coefficient x) then the solution exists globally in time in L^, whereas if the initial mass 
is larger than m* then the solution blows up in a finite time. We refer to the surveys [?, ?] 
for a complete and detailed description of the literature of this topic. 

In the last years, some authors [?, ?, ?] have proposed variants of Keller-Segel type models 
featuring global existence of solutions no matter how large the initial mass is, obtained by 
replacing the linear diffusion term in the equation for the population density by a degenerate 
nonlinear dijfusion term with super-linear growth for large densities. This choice can be 
motivated by taking into account the fact that cells do not interpenetrate (that is, they are full 
bodies with nonzero volume) and therefore diffusion is supposed to inhibit singular aggregation 
effects when the density is very high. We mention here that other authors proposed the use 
of a nonlinear chemotaxis coefficient px{p) which attains the value zero when the population 
density p reaches a fixed maximal value - see for instance [?,?,?]- being this choice motivated 
by the fact that individuals stop aggregating when the density is too high. In both cases, in 
the resultant model overcrowding of cells (concentration to deltas for the cells density p) is 
prevented independently on any initial parameter. 

In the last years, several authors have started to describe biological systems with chemo- 
taxis from a hydrodynamical point of view, i. e. via nonlinear hyperbolic systems of Euler 
type, see in particular [?, ?, ?, ?, ?]. The models obtained are of the form of system (jl.ip . 
where the chemotactical force Vc in (II. ip and the pressure contribute to balance the rate 
of change of the momentum. Moreover, our model (jl.ip features a friction term modeling 
the drag between cells and the substrate material (some authors also considered models with 
a linear viscous term). In this framework, the nonlinear pressure term g{p) in (jl.ip plays 
the role of the diffusive one in the drift-diffusion equation. Therefore, one can interpret the 



overcrowding-preventing effect described before (sometimes referred to as volume filling ef- 
fect) by thinking of the cehular matter as a medium with limited com,pressibility, i. e. closely 
packed cells exhibit a limited amount of resistance to compression. In this sense, a reasonable 
choice of a pressure g{p) is a function of the form g{p) = p"^, ■y > 0. Such an expression also 
has the advantage of modeling absence of stresses for low densities (see [?] for a more detailed 
description). 

In this paper we want to contribute to the problem of establishing a rigorous mathematical 
link between the system (jl.ip and several Keller Segel type models of the form (jl.2p in terms 
of diffusive relaxation limits. A typical example of diffusive scaling on the system (II. ip that 
we shall consider (see subsection 12. ip is the following 

, 1 t ,, , I f t 

a= -, T = -, V [x,t) = -V I X, - 

e e e \ e 

which transforms p.ip into the following rescaled system 

dtp^ + V • {p'^v^) = 

e^ [dtv^ + v^ ■ Vv^] + Vg{p^) = Vc" - v' (1.4) 

^edtc^ = A(f + Q/9^ — j3(f 



r£ 



Formally, as e — > 0, we expect the solution [p^^v^,c^) to system (II. 4|) to behave like the 
solution {p^ ,vP ,c^) to the drift-diffusion system of Keller-Segel type 

ra,pO + v(pOv(cO-5(p°)) = o 

where the loss of the persistence term in the equation for the momentum yields a constitutive 
law for the velocity v^ = Vc^ — 'S/g{p^) (which can be considered as an equivalent of the Darcy 
law in [?]). 

A way to understand the meaning of this phenomena is to consider it as the large time 
behaviour of dissipative nonlinear hyperbolic systems and to look at the asymptotic profile 
as the relaxed equilibrium. This is the case for many relevant situations in mathematical 
physics and applied mathematics. Singular limits with a structure similar to (jl.4p have been 
analyzed by Marcati and Milani [?]. In that paper they investigate the porous media flow 
as the limit of the Euler equation in 1 — D, later generalized by Marcati and Rubino [?] to 
the 2 — D case. Relaxation phenomena of the same nature appear also in the zero relaxation 
limits for the Euler-Poisson model for semiconductors devices and they were investigated by 
Marcati and Natalini [?, ?] in the 1 — D case and by Lattanzio and Marcati [?] in the multi-D 
case. For a general overview of the theory of the singular limits see the survey [?] and the 
paper [?], where the theory is completely set up. 

To perform the relaxation limit we follow the same techniques developed in [?, ?, ?] (among 
others), which are crucially based on the method of compensated compactness of Tartar and 
Murat (see [?, ?, ?]) combined with the Young measures associated to the relaxing sequence 
p'^ (see [?, ?, ?, ?, ?]). Throughout the whole paper, we shall restrict ourselves to the case of 
two space dimensions, which is also the most treated case in the literature concerning Keller- 
Segel type systems. Moreover, for the sake of simplicity we shall work on the 2-dimensional 
torus. We shall prove that this singular limit can be rigorously justified as far as the new time 
variable r stays in a bounded interval [0, T] for an arbitrary T and provided that certain a 



priori assumptions holds for the solution to (jl.4p (see assumption 12.21 below). These a priori 
assumptions are usual in the framework of relaxation limits for nonlinear hyperbolic systems 
(see also [?], [?]) and they don't include any smallness assumption on the initial conditions. 
The rigorous statements of these results are contained in Theorem 14.31 

In order to produce a nontrivial class of solutions to the nonlinear hyperbolic system (II. ip 
which relax toward a Keller-Segel type model after a proper rescaling, we shall also provide an 
existence theorem for the approximating system ()1.4p and prove the uniform estimates needed 
to justify the assumptions (|2.2p in case of initial densities po which are small perturbation 
of an arbitrary non zero constant state (see Theorem 15. 2p . This result is achieved by means 
of the classical FRiedrichs' symmetrization technique and by a linearization argument, see 
[?, ?]. We remark that, in many of the estimates performed here, the pressure term need 
not to be of the form g{p) = p^ . Indeed, some of the estimates proven are still valid if one 
considers a logarithmic pressure g{p) = logp, which corresponds to a linear diffusion term in 
the limit problem ()1.2p (this fact is not in contradiction with the blow-up of the density in 
the limit problem with linear diffusion, see the Remark 15. 4p . However, while considering the 
alternative scaling introduced in subsection 12.21 (where the limit is the classical Keller-Segel 
system (|1.3p ). such an expression for the pressure seems to be essential in order to achieve 
the needed estimates no matter how large the initial mass is, in a very similar fashion to what 
happens in [?]. We remark that our convergence results hold on an arbitrary time interval. 
Therefore, at least in the case of the second scaling treated in 12.21 (where the expression 
g{p) = p'^ is crucial in order to achieve the desired estimates), our result can be seen as a new 
interpretation of the overcrowding-preventing effect due to the power-like expression of the 
pressure. More precisely, the global smoothness of the limit density p^ (and the absence of 
concentration to deltas for all times of p^ as a byproduct) can be obtained as a consequence 
of our relaxation result, alternatively to the more direct proof developed in [?, ?]. 

The paper is organized as follows. In chapter [5] we state the three different scalings we 
shall deal with. In chapter [3] we perform the uniform estimate we need in order to prove the 
main convergence theorem. In chapter U] we prove the main convergence theorem for large 
data under the a priori assumption 12.21 by means of compensated compactness and Minty's 
argument. In chapter [5] we prove an existence theorem for the approximating rescaled system 
in order to provide a class of solutions satisfying the basic assumptions ()2.2p . 

2 Preliminaries and rescalings 

We consider the following •persistence and chemotaxis model 

drP + V- {pv) = 

drV + v-Vv + Vg{p)=Vc-dv (2.1) 

adrC = Ac + ap — j3c 

where a, /3, d, a are nonnegative constants. The system (|2.ip is endowed with the following 
1-periodic initial data 

p(0,x) = /)o(x), v{Q,x)=vq{x), c(0,x)=co(x) 

The nonlinear function g[p) has the form 

g{p) = p'^ , for some 7 > 0. 



Remark 2.1. The nonlinear function g{p) grows faster than K*\ogp, for large p, where 
K* = M/Att and M is the total mass of p. More precisely, there exists U > 0, such that 



for any p>hl and k > k* g{p) > nlogp, 



(see [?]). 



Some of the results contained in the present paper hold in any space dimension n, whereas 
some of them are true only in the case n = 2. In order to simplify the coverage, we shall 
always restrict ourselves to the latter case. In the sections [21 E] and S] we shall not deal with 
the existence theory of ()2.ip . whereas we shall work under the following basic assumption. 



Assumption 2.2. There exists a global solution (p, v, c) to ()2.ip . smooth enough in order to 
justify the estimates contained in section [3] and such that 

(Al) the total mass M = J pdx is conserved, 

(A2) p{x,t) >k>0 

(A3) {p,pv)eL°°{T^x [0,+oo)). 

Let us now explain in detail the relaxation limits we want to perform. We shall deal with 
three different asymptotic regimes for p.ip . corresponding to small parameter limits of three 
different types of scaling. 

2.1 First scaling: large times and large damping 

For a fixed constant e > we consider the large damping rate d = - in (12. ip . Then, we 
introduce the fast time variable 



t 

e 



and the new independent variables 



f^(x, t) = -t;(j;, r), p'^{x,t) = p{x,t), c^{x,t) = c{x,t). 



(2.2) 



Moreover, we fix cr = 1 in the third equation. Then, system (12. ip in the new variables reads 

(2.3) 



dtp" + V • [p^v") = 
£2 [dtv" + v^ ■ Vv"] + Vgip") = Vcf 
^edtc" = Ac" + ap" - (3c". 



The formal limit as e ^ is given by the parabolic-elliptic system 

fp? + V-(p0v(c0-5(/50)) = 
1 Ac° + ap^ - I3c° = 0. 



(2.4) 



2.2 Second scaling: large time and large damping with Poisson coupling 

A simplified version of (12. ip . namely with (3 = and a = 0, is given by the following system 



drP + V- (pv) = 
drV + V ■ Vv + Vg{p) = Vc — jv 
^0 = Ac + ap. 



By performing the same scaling as before, namely 



r = -, v'^{x,t) = -v{x,t), p'^{x,t) = p{x,t), d^{x,t) = c{x,t), 



and by putting d = j, we obtain 



'dtp^ + V ■ (p^t;^) = 
= Ac^ + a/)^ 



(2.5) 



(2.6) 



(2.7) 



The formal limit as e ^ leads to the usual Keller-Segel model with nonlinear diffusion 
(see [?]) 

\ AcO + apO = 0. 



(21 



2.3 Third scaling: diffusive scaling with small reaction rates 

Starting once again by ()2.ip . we consider the case a = d = 1. We consider e-depending 
reaction coefficients a and (3, namely we require a = ea and /? = e/? for fixed 5, /3 > 0. We 
then perform the diffusive scaling 

t X 1 

T = — , y = -, v^{x,t) = -v{x,t), p''{x,t) = p{x,t), c^(x,t) = c(x,r). (2.9) 

This leads to the rescaled system 



dtp" + V • (p^i;^) = 
e^ [dtv^ + v^ ■ Vv"] + Vc/(p^) = Vc" - v" 
^ dtc" = Ac" + ap" - Pc. 



(2.10) 



Therefore, the formal limit as e — > is given in this case by the following fully parabolic 
model (we drop the ~ symbol for simplicity) 



pO + V-{p^V{c^-g{p^)) = 
qO = AcO + apO - /3cO. 



(2.11) 



Remark 2.3. From the hypotheses (A3) and the scalings (f2^ . (f2l6]) . (fZ^Ql) . it follows that 
the sequences {p"}, {ep"v"} are uniformly bounded in L°°(T^ x [0,+oo)) with respect to e. 



3 Estimates 



In this section we provide suitable estimates on the solutions of the three rescaled models 
([23]) . (I2TII and (l2:T(l . For future use we define 



rp I 

Pip) ■= / 9{n)dn = — — ; 
Jo 7+1 



,7+1 



(3.1) 



3.1 First scaling 

We have the following (standard) energy estimate for the rescaled system (12.30 . 

Proposition 3.1. The following identity is satisfied for any t G [0,T], by any solution 
(p^^;^c^) to fH 



-p'{x,t)\v'ix,t)\^+P{p%x,t)) 



dx -\ — / / p^ {x , s)\v^ (x , s)\'^ dxds 
2 Jo Jt^ 



-p§(x)bg(x)p + P(pg(x)) 



dx + [ Kt + e 



c5(x) 



dx 



J2 



(3.2) 



Proof. By multiplying second equation in (|2.3|) by p^v^ by using the first equation in (|2.3|) 
and by integration by parts it follows that 



d_ 

dt .ij2 



£ 



■p'{x,t)\v'{x,t)\''+P{p'{x,t)) 



dx + p'^{x,s)\v'^{x,s)\dx = 

Jt2 

p'^{x,s)v^{x,s)Vc'^{x,t)dx < - / p'^{x,s)\v^{x,s)\'^dx + -\\p^\\oo / \V c"^ {x , t)f dx . 

J2 2 Jj2 2 Jj2 

(3.3) 
Now, by multiplying the third equation of (j2.3|) by c^ we get for any 6 > 

u / fc, -, .,,0 , / , — ,- , .,,0 , ,, ^,, / il 



dt 



|c^(x,t)rdx 



■ / |Vc^(x,t)p(ix + a||/)^ 

Jt2 



00 I -J^+^ l\c%x,t)\^dx 



P / \c''{x,t)\^dx 

Jt2 



(3.4) 



By choosing 5 < f' by integrating in time we obtain, for fixed constant K, independent on 
e, that c^ satisfies the following inequality 



e. - .m2j^ , /5 



c%x,t)\''dx + '- I [ \c%x,t)\'^dxds+ [ [ \Vc%x,t)\^dxds<Kt + e [ *^°;^* dx. 
T2 2 2 Jo JT2 Jo Jt2 Jt2 2 

(3.5) 

The estimate p.2p follows now by using together ()3.3p with ()3.5p and by taking into account 

the hypothesis (A3). D 

Corollary 3.2. Lei {p^,v^,d^) be a solution to 112. 3\) satisfying assuniption \2. 2\ and with initial 
datum {pq,Vq,Cq) satisfying 



-pl\v^^\^dx + {pf^r+' + e\cl\' 



dx uniformly bounded w.r.t. e <^ 1. (3.6) 



Then, for all T > 0, 

Ey/^v^ is uniformly bounded in L°°{[0,T],L^{T )), for all p > 1, (3-7) 

p' is uniformly bounded in L°°([0, T], ^^(T^)), for all p>l, (3.8) 

^fffv'^ is uniformly bounded in L^([0, T] x T^), (3-9) 

y/ec" is uniformly bounded in L°°([0, T],L'^{T'^)), (3.10) 

c^ is uniformly bounded in L'^{[0,T],H\t'^)). (3.11) 

Proof. ()3.7p and ()3.8p are a consequence of the assumption 12.21 while (13. 9|) follows from the 
inequahty ([32]) • Finally (fXTO]) . (f3lT]) follow from (l33|). D 

3.2 Second scaling 

We consider the following energy for the solution to ()2.7p 

^.(0= / Kp'\v'\''+p{pn-lp'c' 

./in>2 z 2 



dx, (3.12) 



where P is given by (|3.ip . For semplicity we will take a = 1. We have the following estimate. 



Proposition 3.3. The following inequality is valid for a solution (yO^,u^,c^) to (2.1): 

Ee{t)+ / p^x,s)\v^x,s)\'^dxds = Es{0) (3.13) 

Jo Jt2 

Proof. By using the Poisson equation of the system ()2.7p we easily have 
— ^^(t) = - / p^|t;^|2dx+ / p^t;^-Vc^dx+ / c'Acfdx 

at Jf2 Jj2 Jj2 

= - I p'\v'\^dx+ I p'v' ■Vc''dx+ I {c'V ■{p'v'))dx, 

and this implies 

^E,{t) = - [ p^\v'\^dx. 
at Jj2 

Integration with respect to time yields (|3.13p . D 

In order to recover an estimate for Vc*", let us introduce the following convex functional 

J[P']= [ {P{p'{t))-p'c%t))dx. 



JT2 

Now we proceed by estimating the functional J[p^] from below, using the same strategy of 
[?]. Let us recall that if c^ S VF^'"'^(T^), then the convex functional J[p^] has a critical point 
p* which is a solution of 

g{p*)-c' = X (3.14) 

whenever p* > and null otherwise. Here A is the Lagrange multiplier associated to the 
constraint given by the mass conservation J p* = M and fixed by this condition. We refer to 
[?] and ([?], Proposition 5) for details. Therefore, we have 

J[p'] > [ {P{p*) - P*c')dx = [ {P{p*) - p*g{p*) + \p*)dx. 

JT2 J{p*>0} 



By taking into account the Remark 12.11 we can introduce the corrective term R such that 
g{p*) = Klogp* + R{p*), then we have 

J[p']> I iP{p*)-Kp* log p*)dx- f p*R{p*)dx + \M. (3.15) 

Jf^ i{p*>0} 

Now, p.l4|) imphes nlogp* + R{p*) = X + c^ , whenever p* > and thus 

[ exp f :^^^ p*dx = e^/" / exp ( -) dx, 

^{p*>0} \ '^ J J{p">0} \i^J 

so we have 

A = Klog ( / e^^''p*dx] -Klogi [ e^'/'^dx] . (3.16) 

\J{p'>0} J \J{P'>0} J 

If we replace A by its expression in the inequahty ()3.15p , we conclude that 
J[p']> f {P{p*)- up* log p*)dx- [ p*R{p*)dx 

+ kM log I [ e^^''p*dx] -kM log I [ e^'/^dx] . (3.17) 

\J{p*>0} J \J{P*>0} J 

By taking into account the Remark 12.11 we have that 

f {P{p*) - Kp* log p*)dx > C. 

J\p*>U\ 



l{p*>U} 
On the other hand, we have 



/ {P{p*)-Kp* log p*)dx> - ( sup(P-Kplogp)" I |t2|. 

J{p*<U} \[a,U) J 

Therefore, 

/ {P{p*)-p* log p*)dx 

is uniformly bounded form below. Now, the Jensen inequality for the probability density 
p* jM over the set where p* > 0, gives us that 



{p*>o} \ K M J J{p*>o} M 



exp I { ^dx ) < / e^i'^^dx, 

ans thus 



KMlog ( / e^/'^^dx I - / p*R{p*)dx > 0. 

\J{p*>o} M J 7{p.>o} 

Finally, we recall and use the Trudinger - Moser inequality [?, ?, ?]. 

Theorem 3.4. Assume that fi C M^ is a C^, bounded, connected domain. It exists a constant 
Cq, such that for all h G H^ with J^h = we have 

/ exp{\h\)dx < Cuexpl— / \Vh\^dx] . 



By applying the previous theorem to (f j k we obtain 

I e'^'/'^dx < I e^'/^dx < exp ( -^ I \Vc^\'^dx 



'{p*>0} 

and thus 



'{p*>0} 

So by (|3.17|) we have that 



KMlogi I e^'I'^dx] > — ^ [ iVc'Pdx. 



J[p']>C-^ I \Vc'fdx (3.18) 



•2 



Proposition 3.5. Assume {p^,v^,d^) be a solution to 1^2. 7\ ) satisfying assuniption \2.2\ then 

I IVc-r^d. ., umfor^ly .„„„« (3.19) 

Proof. We can rewrite p.lSp as 
^^(0)= / t^p^\v^\^dx+ f J[p'] + \ [ \Vc'{t)\^dx+ f [ p%x,s)\v'{x,s)\'^dxds. (3.20) 

JT2 2 Jj2 2 Jj2 Jo JT2 



Combining I^3l0\i with (l3l^ we get that 

-2 rt 



Ee{0)> %-p''\v^\^dx+ i / p''{x,s)\v''{x,s)\^dxds 



Finally, Remark 12.11 implies k > k*, i.e. (l — ^r) > and thus 

is uniformly bounded. D 

Corollary 3.6. Let {p^,v^,d^) be a solution to \2.T^ satisfying assuniption \2. 2\ and with initial 
datum {pq,Vq,Cq) satisfying 



'lp^^\v^^\'dx+{pt,r+'-lpf,4 



dx uniformly bounded w.r.t. e <C 1. (3.22) 



/T2 

Then, for all T > 0, 

e^v'' is uniformly bounded in L°^([0, T], ^^(T^)), for all p>l, (3.23) 

p^ is uniformly bounded in L°°([0, T],LP{T^)), for allp>l, (3.24) 

\f^v^ is uniformly bounded in L ([0, T] x T ), (3.25) 

c= is uniformly bounded in L'^{[0,T],H^{T'^)). (3.26) 

Proof. ()3.26p follows from Proposition 13.51 and by taking into account that we are in a pe- 
riodic domain. (j3.23p and (j3.24p are a consequence of the assumption 12. 2 1 while (j3.25p is a 
consequence of (I3.13P and ()3.19p . D 

10 



3.3 Third scaling 

With the same procedure as in the Proposition 13. II we are able to prove that 

Proposition 3.7. Let {p^,v^,(f) he a solution to 112. 10\) satisfying assumption \2.S\ and with 
initial datum {pq,Vq,Cq) satisfying 



/ 



pl\vl\^dx + {plV^^ + \cl\'' 



dx uniformly hounded w.r.t. e <C 1. (3.27) 



Then, for all T > 0, 

ey/Jfv^ is uniformly bounded in -L°°([0, T],LP{T^)), for all p>l, (3.28) 

p^ is uniformly bounded in L°°([0, T],LP{T^)), for all p>l, (3.29) 

\ff^v^ is uniformly bounded in L ([0, T] x T ), (3.30) 

c^ is uniformly hounded in L'^{[0,T], L'^{T^)) n L'^{[0,T], H\T^)). (3.31) 

4 Strong convergence 

This section is devoted to the study of the relaxation of the systems (12. 3p , (12. 7p , (I2.10p towards 
their formal limit (j2.4p , (j2.8p , ()2.1ip , respectively. As a consequence of the Corollary 13.21 and 
the Propositions 13.51 13.71 we have that, extracting if necessary a subsequence, 

Vc^ ^ Vc° as e i weakly in L^{[0, T] x T^) . 

This convergence for c^ is enough to pass into the limit in p.3p . ()2.7p . ()2.10p to get in the 
sense of distribution (j2.4p , (|2.8p , (j2.1ip , respectively, p rovided that p^ converges in a strong 
topology. In fact by the Remark 12.31 we know that p'^ — > p^ in L°° *— weakly, while by ()3.8p . 
()3.24p . (j3.29p we have p^ -^ p^ weakly in L^, for any p > 1. These convergence are clearly too 
weak to pass into the limit in the nonlinear terms of ()2.3p . p.7p . ()2.10p . So, in this section 
we will investigate the strong convergence of the approximating sequence p^. The analysis of 
this convergence reduces to the analysis of the convergence of quadratic forms with constant 
coefficients via the classical compensated compactness technique due to Tartar [?, ?] and 
Murat [?] (see Dacorogna [?]). As we will see later on, these techniques will apply in the 
same way to the three scalings (|2.2p . (j2.6p . (j2.9p . so we will discuss them together. Let us 
recall the following theorem 

Theorem 4.1. (Tartar's Compensated compactness) 

Let us consider 

1. a bounded open set Vt C M"; 

2. a sequence {l'']f=i, l" : QcW — > MT; 

3. a symmetric matrix G : M™ — > W^ ; 

4- constants a^,-^ G M, i = 1, . . . , g, j = 1, . . . , ?n, A; = 1, . . . , n. 
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Let us define 



f{a) = (0a, a) , for all a G 
A= < 



A G R™ : 377 G R" \ {0}, J^ aj.A,^ = 0, i = 1, 



> . 



Assume that 



(a) i/iere exists I G L^ (il) suc/i t/iai l'^ ^ I in L^ (i7) as 1/ 1 00; 

(b) A^l'^ = Yl ^\k'dx~' ^ = 1) • • • <? OjTS- relatively compact in H^^^ (Q); 

(c) /|A = 0; 

(d) there exists / G R such that f{l) ^- f in the sense of measures M{0,). 
Then we have f = f{l). 

4.1 Weak convergence of p'P^p'^) 

First of all we start by studying the weak convergence of p'^P{p'^). Our goal will be to prove 
that 

p'P{p')^p'P{p'), 

where p" is the weak limit of p"^ . To this end we are going to apply the Theorem 14.11 in the 
same spirit of [?]. In order to fit the into the hypotheses of the Theorem 14.11 we rewrite the 
first two equations of the systems (|2.3|) , (|2.7|) , (J2.1U|) , as 



pf + m| + n^ = 



e ml + \ e 



e^t + i E^ 



rWu' 



P' 



+ lP{p')\ +{e 



P' 



[n 



s\2 



+ e'^^ + ^P{p' 



P'cl - Tjf 



p'^cl - n^. 



(4.1) 



P J X \ P" 

where 

v^ = (wf , vf) p^'v" = ip^vi,p^V2) = {m^,n^). 

It will be usefull rewrite (14. ip in the following way 



(4.2) 



lPip% 



-e^ml - e^ 



,m 



e\2 



P' 



e 



^P{p% = -e-n\ -e-\'^ "^ 



p- 

.2 / {n^' 



+ p cj. — m 



(4.3) 



X ^ ' ■' y 



P^ 



+ p'ci - n^ 



y 



By using i^M), dSSS]), ([3301), ([XT^ . and the assumption (A3) we get that p^v^ p^Vc^ G 
L2([0,r] X T2). In fact 

ll/5^^'''llL2([0,r]xT2) < ||\/p^||oo||\/p^v''||L2([o,T]xT2) (4-4) 

IIP^Vc^IIl2([0,T]xT2) < ll/'^l|oo||Vc^||i2([o_T]xT2) (4.5) 
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Moreover, by taking into account the assumptions (A2) and (A3) and (j3.9p . (j3.25p . (j3.30p we 
have that e^ I ^™ ^ I , is relatively compact in H^^{[0,T] x T^). In fact let us consider co 

relatively compact in [0, T] x T^, then by taking into account (A2), (A3) and the Remark 
we have, 



[m 



e\2 



H-iH 



< sup 



(m^) 



e\2 



1 



< e||/9^7;'||oo^||\/p^w1lL2(^) 



2 / m'^n'^ 



p- 



, e 



2 / va'^n' 



In a similar way it can be proved that e^ 

relatively compact in H-^{[Q,T] x T^). Now, (I13D-(jMD, imply that 



2 ({n-f 



'Pt +m% + nl 
P{p% 
P{P% 



n^ 



is relatively compact in (H^^^) 



(4.6) 
e^{p'^v'^)t are 



(4.7) 



In order to fit into the framework of the Theorem 14.11 we set xi = x, X2 = y, x^ = t, 
F = {m^,n^, p^,P{p^)), hence m = A. In our case the differential constraints are q = 2>. So 
we can define the matrices A^,J?,A^ G A^4x3, where A^ = {a*fc}, ^ = !> 2, 3, j = 1, . . . ,4, 
A; = 1, 2, 3 as follows: 



A^ 



(I 0\ 
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Vo 0/ 



A^ 



/O 0\ 





\l o/ 



A' 



'jki 


1 <- 


— 1 


/o 





0\ 




















Vo 


1 


0/ 



The characteristic manifold A is then given by 

A = {A G M^ I 3e G M^ \ {0}, Bit A) = O} 



where 



Therefore 



If we define 



'AiCi + A26 + AsCa 
B{tX)=\ A4ei 

A4^2- 



A= <^ 



Ai A2 A3 ' 
A G M"^ I det I A4 
A4 



{AgM^ I A3A4 = 0} 



2 



/O 0\ 



£M 



4x4) 




1 

Vo 1 0/ 

then /(A) = A^MA = A3A4 and, of course /|a = 0. Now, by applying the Theorem 14.11 we 
have f{l^) -^ f{l), and in our case this means 

p'Pip') ^ p'P', 
where P^ = w — lim.P{p^). 
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4.2 Strong convergence of p^ 

In the previous section we proved that 

Here we wih be able to prove that 

p^ -^ p° strongly in Lf^^, p < +00. 

At this point we can follow the methods of [?], [?]. First of all let as use Minty's argument 
([?], [?]) to prove that P^ = P{p^). Since the function P is monotone, for any w G L°° e (f 
test function, 99 > 0, we have that 



H{e) = {P{p') - P{w)){p' - w)^dxdt > 0; (4.8) 

But for e I 0, we have that 

f f Pip')p''ipdxdt^ f f P^p^^dxdt, 
So from ()4.8|) we get that for e | 

H{e) ^ H= f f{P^ - P{w)){p^ - w)^dxdt > 
If we choose w = p^ + Xz, with A < and arbitrary z G L°°, we have 
G(A, z) = f f{P° - P{p^ + \z))z^dxdt 

= i [ [iP°- PiP° + ^z))Xz(pdxdt < 
and for A T 0, G(0, z) < for any z G L°°, then 

^(0, z)= f f{P^ - P{p^))z^dxdt = 0, 



and finally P° = P{p^). 

Our next step now, is to prove the strong convergence for p^ — > p^ in L^^^. To this end we 
characterize the weak convergence by means of Young's probability measures (see [?] , [?] , [?] , 
[?], [?]). Let us recall that if {u^} is sequence converging to U in L°° *- weakly , we can 
associate to {u^} a family {i^(x,t)W} of probability measures such that for any continuos 
function F{-) 

* - Imi F{u') = J F(A)z.(,,i) (dX) a.e 

If ^(x,t) = ^U{x,t) then u'^ ^ U strongly in L^^^ for any p G (l,+oo) (see [?], Corollary 6.2). 
Let {i^(x,t)} be the family of Young's probability measures associated to the sequence {/O^}: 
since p*" — > p^ in L°° *- weakly, we can find a closed interval [a,b], < a < b, such that 
supp i'ix,t) — [o, ^]- Since P{r) = r", a > 1, we have three possibilities: 
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1. P G C2 (M \ {0}) e P"ir) j +00 for r i 0, if 1< a < 2; 

2. P £C^ (E) e P"(0) = 1, if Q = 2; 

3. P £C^ (R) and P"(0) = 0, if a > 2. 

Let us assume that 1 < a < 2. Then we can write for any A, Aq 

P(A) - P(Ao) = P'(Ao) (A - Ao) + ^P" (A*) (A - Aq)' , 
where A* belongs to the segment between A and Aq. If we choose 

Ao = / Az/(^^i) (dX) = po, 



since P° = P(p°) 



so that 



P(Ao)= / P(A)zv(,,t)(dA) 

J a 
b 

{P(A) - P(Ao)}i^(.,i) (dX) = 0. 
On the other hand we also have 

P'(Ao) (A - Ao) z^(,,i) (dX) = P'(Ao) 1^ Xu^^^t) (dX) - Aq ^ u^,,t) idX)\ = 0, 

so we can conclude that 

' P"{X*){X-Xofu^.,,t){dX) = 0. 

Taking E = min P" (A) > 0, we get 

Xe[a,b] 

f-b 

E {X- Ao)' z/(,,i) {dX) < 0, 

J a 

namely 

(A-Ao)2z.(,,)(dA)=0, 

and it follows a = b and vix^t) = (5) a point mass and so we finally get 

p^ -^ Po strongly in Lf^^. 

To conclude we remark that in the case a > 2 this result can be obtained in the same way 
by using the function — P^^. 

Remark 4.2. The strong convergence result for p^ obtained in this section is still valid in the 
case of linear diffusion, namely if we consider g{p) = logp and consequently P{p) = plogp— p. 

By using the estimates and the strong convergence of the sequence {p^} obtained in the 
previous section we are able to prove the following main theorem. 
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Theorem 4.3. Let T > be arbitrary and let {p^,v^,c^) be a family of solutions to the 
system (I23D TdSI) and ([TTU]) respectively) with initial data satisfying (EB* TlEU) and [^^FJ\ ) 
respectively). Assume that the assumption \2.^ holds, then, there exist p^ G L°°([0, T] x T^) 
and c^ G L^([0,T],i/-'^(T^)), such that, extracting if necessary a subsequence, 

/9^ — > p^ strongly in L^^^([0, T] x T^) for any p < oo 
Vc" -^ Vc° weakly in L^{[0,T] x T^). 

Moreover the couple {p^,(P), satisfies the system ()2.4|) (^ (I2.8p and ()2.1ip respectively) in the 
sense of distributions. 

5 Perturbation of constant states in the approximating system 

In this section we deal with the rescaled system 

dtp" + V • {p"v") = 

£^ [dtv" + v" ■ Vv"] + Vg{p") = Vc' - v" (5.1) 

edfc'^ = Ac^ + ap^ — /3c^. 

with X G T^, t > 0, where T^ is the flat normalized two-dimensional torus. The system (j5.ip 
is complemented with the 1-periodical initial data 

p'{x,0) = pI{x), v'{x,0) = v'oix), c'{x,0) = clix). 

We shall consider small perturbations of the stationary state 

^ ^ ^ ^ ^ ex ^ 

{p,v,c) = {p,v,c}, p>0, v = 0, c=-/9 (5.2) 

and prove the existence of solutions {p^,v^,c^) such that the density p^ stays away from 
zero, uniformly in e, on a small enough time interval [0, T] with T independent on e (see 
similar results in [?] and [?]). In order to perform this task, we shall use an iterative method, 
namely we define recursively the sequence (p", v", c") as follows: {p^{x, t), v^{x, t), c^{x, t)) = 
{Pq{x),Vq{x),Cq{x)) and, for n > 1, (p"',i;",c") solves the linear system 



dtp"" + «"-^ • V/5" + p^-^V • n" = 

,n-l 



0'f/9"~l) 1 1 



dtc"" = -Ac" + -/>" - -c". 

£ £ e 

From now on we shall drop the dependency on e on the solutions (p, v,c) to simplify the 
notation. Moreover, we shall use the following notation: the variables taken at the step n — 1 
will be denoted e. g. by p^~^ = p; the variables taken at the step n will be denoted without 
any further symbol, e. g. p" = p; the deviation from the aforementioned constant stationary 
states will be denoted e. g. by p = p" — p and p = p"^^ — p. 

The first two equations in system ()5.3p can be easily viewed as a hyperbolic system in 
vectorial form. More precisely, let us define the 3-dimensional variable U as 

U := ip,v\v'), 
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where v = {v ,v ). Let us denote 



A,{U) :-- 



{v' 


p 


01 


g'ip) 


v' 











^^J 



A2{U) :-- 



v-" p 
7)2 

^ ^2 



B{U) :-- 



Then, with all these notations, the system (15. 3p can be rephrased as 



dtU + +Ai{U)d,,U + A2{U)d,,U = B{U) 

Otc = -Ac H — p c. 

£ e e 






(5.4) 



The first line in (|5.4|) corresponds to a linear hyperbolic system which can be symmetrized by 
means of the matrix 



5(C/):=diag(^,p,p 



(5.5) 



The matrix S{U) is uniformly positive definite provided the variable /) satisfies a condition of 
the form < c < p < C (we recall that g'{p) = 'yp'^~^ exhibits a singularity at zero in case of 
7 < 1). The two matrices S{U)Ai{U) and S{U)A2{U) can be easily proven to be symmetric. 
We now rewrite system ()5.4p in terms of the deviations U and c: 



dttJ + Ai{U + U)d^, U + A2{U + U)d^,U = B{U + U) = B{U) 

OfC = -Ac H — p c. 

e e e 

We introduce the energy functional 



(5.6) 



8{U,c) 



1 



U' S{U)U + \c' 



dx 



1 



g'ip) 



2 I ^1 |2 I \ 2 

— p +p|t;| + Ac 



dx, 



where A > is a constant to be chosen later on. We have the following 

Proposition 5.1. Let T > 0. There exist constants eo,6 G (0,1), K G {0,p/2) such that, if 

WPo - pIIh4(t2) + e||^ollH4(T2) + \/e||co - c||^4(Tr2) < 6 and 



sup {\\p{t)\\Hi(j2) + e||i;(t)||//4(T2) + \/e||c(t) - c\\Hi(^f2)) < K, 

0<t<T - '- ' 



for all e G (0,eo); then, 



sup 

0<t<T 



H4(T2) + e\\v{t)\\H'i(T2) + Ve\\cit) - c||H4(Tr2)) < K 



(5.7) 



(5.8) 



for all e G (0,eo) 



Proof. During the proof of this proposition we shall often make use of the Sobolev inequality 
||/||Loo(Tr2) < C||/||j^2cir2). 

Step 1. Due to the symmetry of the two matrices SAi and 5^2) we can use integration 
by parts in the evolution of £{U, c) as follows: 



l^'^'^'^i- 



T 



SiU)AiiU) 



1 



Udx + - U 

XI 2 Jj2 
T otTT\T:>ti~T\J ^ f IV7-|2 7 , "^^ f 7 ^(^ 



tT 



S{U)A2{U) 



Udx 



X2 



+ / U' S{U)B{U)dx - - |Vc|^dx + — / pcdx - -- 

Jt2 £ Jj2 e Jf2 e Jj2 



c dx. 



(5.9) 
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Due to the assumption (|5.7p we have 

7 ^ P y -2 \—\2\ \\^\ r 

-£{U, c) < C{K) (IIVpIIloo + \\Vv\\lo^) 2 / J ^ + T^j ^^ + ^^^ / / • ^^^^ 

_ (p-_^ /■ |^|2rf^_^ / |Vc|2dx + — /" /^cdx-^ /" c?dx (5.10) 

for a function C(K) > of the constant K such that C is continuous on ii' G [0, /)/2]. By 
choosing 

_. JP + Kf 
e{p-K) 

and eo < 1) we can use once again (|5.7p and find a constant Ci > such that 



|.,.. , . .c,A-,i /^ (g . ^) .. - <^ 1^ Kt.. - J|±|j^ X^ ,v.r^^ 



(p + K)2/3 /■ _2^ , (p + K)V / _2 



2(p - i^)e 7,^2 ' '^'^ ^ 2(p - K)e/3 7^^. ^ ^^^ 

We now choose the constant K such that KC{K) < 2^ ~ ^)- ^y using the coercivity 
property 



£{U,c)>c{K) f 

7X2 



A' I i9 ^ 

^ + \v\^ + — 
e^ e 



dx, (5.11) 



which holds for a certain c{K) > 0, due to Gronwall inequahty we easily obtain 
£(U(t),c(t))+\ f I nu(r)|2dx + |Vc(r)Pld2;dr+- / I c{TfdxdT<A£{UOd),cme^^ 

£^ Jo 7x2 e 7o 7x2 

for certain constants A,B > depending only on K and Eq. The above implies in particular 
sup / [pitf + e>(t)P + ec(t)2] dxdt < C{K)6e^^, 

0<t<T Jt2 

for a certain C{K) depending on K. Therefore, by choosing 6 small enough such that 
C{K)5e^'^ < K"^ we obtain 

sup (||p(t)||L2(Tr2) +e||t;(t)||L2(Tr2) + Ve||c(t)||L2(T2)) <K. 

0<t<T 

Step 2. We now perform the energy estimate of the space derivatives of {U,c). For 
j = 1, 2 we denote the derivatives with respect to Xj by the subscript pj = dx^p- The system 
satisfied by the first derivatives of (C7, c) is 

'dtUj + Ai{U + U)dx,Uj + A2{U + U)dx,Uj = B{U)j - Ai{U + U)jUi - A2{U + U)jU2 

dtCj = -Acj + jpj - -Cj. 

(5.12) 
The evaluation of the energy 



£iUj,Cj) = \ I \uJS{U)U,+Xc] 



dx 
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in a similar way as in (j5.9p yields 



^ -'- - ^ ^ ^ -l-T [c/m /i^ml T-T.j^ , 1 / t'tT 



£{U„c,) = - U SiU)A,iU) U,dx + - U SiU)A2iU) 

dt Z J'f2 ■> I J Xi Z Jf2 ■> I 



X2 



Ujdx 



TTanJ^T,fT-T^ j^ ^ f iw. |2^^ , ^" f . . ^. ^/^ / .2, 



+ / Ul S(U)B(U)idx / |Vc,rdxH / pjCidx -/ c^dx 

UfS{U)Ai{U)jUidx- I UjS{U)A2{U)jtJ2dx. 

Assumption (j5.7p allows for the estimate of the first two terms above as in (j5.9p . as well as 
for the estimate of the last two terms in a similar fashion. The result is the following estimate 



d 



-rS{U,rc,) < C{K) {\\Vp\\l^ + \\Vv\\l^) - -| + ^ dx + 



dt 



P? 1^,-P 



2 J'f2 \ e^ 



(p-K) 



A 



Xa 



Vj\ dx / |Vcj| dx H / pjCjdx / cdx 



£ Jj2 



\P\\l° 



xp 



Vj ■ Vcjdx 



e Jj2 



e Jj2 



which is the equivalent of the estimate (IS.lOp where {U,c) are replaced by their first deriva- 
tives. Therefore, we can easily conclude as before 



sup (||V/j(t)||i2(Tr2)+e||Vt;(t)||i2(Tr2) + VF||Vc(t)||i2(Tr2)) < K. 

0<t<T 

Step 3. The second space derivatives of {U,c) satisfy the system 

dtUij + Ai{U + U)d,,Uij + A2{U + U)d,,Uij = B{U)i^-Ai{U + U)iUi^-A2iU + U)iU2j 

-Ai{U + U)i,Ui-A2{U + U)ijU2 
-Ai {U + U)jUij -A2{U + U)jU2j 

r.. _!._ ^a_ _P_ 

(^t^ij — ^Cij -]- Pij ^ij 1 

(5.13) 
for i,j = 1,2. The structure of system (I5.13P is similar to (I5.12|) and therefore the estimate 
of the energy 



v^u ' ^ 



ij , ^i] j 



J2 



Ul,S{U)U,j + Ac! 



IJ ^ '^^l] 



dx 



can be performed as in step 2. The only extra terms which needs to be analyzed are the 
following, for i,j. A; = 1,2 [C{K) denotes a generic constant depending on K): 



Uf^S{U)Ak{U),jUk < ^^ / mj\ [\UiUj\ + \Uij\] \Uk\dx 

j2 e Jj2 



< 



< 



K^C{K) 



\Uij\\Uk\dx + ^^^ [ \U,j\\UiAdx 

j2 e^ Jj2 



{K + K^)C{K) 



\U^i\^dx + K^ 



where we have used once again (j5.7p and the result in step 2. Notice that so far we have used 
L°° estimates only up to the firs order derivatives of U and U_. In the last inequality above, 
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the second derivatives are only estimated in L . We have therefore obtained, for Q < K < 1, 






£2 



/" , ,2 , A /• ,^ ,2 , Aa /■ _ , A/3 /• 2 , 
/ l^jjl ax / IVcjjl ax H / pijCijdx / c^jax 



and, by using the same choice of A and K as in step 1, after using Gronwall Lemma we obtain 

£{Uit),c{t))<C{K) 



£:(c7(o),c-(o)) + ^ 



e^*. 



Then, the coercivity property (jS-llh and the assumptions (15.70 imply 

sup iWD^mWLHT^) + e||I?'t^(t)||L2(T2) + V^WD^mh^T^)) < C7(i^)(<5 + K^) 

0<t<T 

and clearly, a choice of 6 and AT small enough implies C{K){6 + K^) < K^, which concludes 
the estimate of the second derivatives. 

Step 4. In order to conclude the proof of the proposition, one needs to perform the 
same energy estimate also on the space derivatives of order 3 and 4. All the estimates on 
the nonlinear terms on the right-hand side are analogous to those in Step 3. The integrals 
with over-quadratic terms always contains not more than two terms involving more than two 
derivatives. Therefore, all the extra terms can be estimated in L°° by using assumption (15. 7p 
and the results in the previous steps. We shall skip the details of these computations. The 
proof is complete. D 

We are now ready to state the main theorem of this section. 

Theorem 5.2. Let T > and let < s < 4:. Let [p,v,7^ he the constant state in \5.S^) . There 
exists constants (5, eo £ (0)1) such that, if the initial data pQ,vo,co satisfy 

WPo - p\\h^{j2) + e||wollH4(T2) + y/e\\cQ - 011^4(1-2) < 6, 

for all £ G (0,eo), then there exists a classical solution (p^,v^,c^) to ( fj.ij) such that the 
quantity 

sup (||p^(i)lk-(T2) + e\\v%t)\\Hs(j2) + ^/e\\c''{t)\\Hs(^f2)) 

0<t<T 

is uniformly bounded with respect to e £ (0,eo) and such that the density p^ satisfies 

p'{x,t) >p/2>0 

for all £ G (0,eo)- 

Proof. For any fixed £ G (0, eo), the sequence (p", u", c") has all space derivatives up to order 
4 in L^ and all time derivatives up to order 3 in L^. Therefore, (/o"',^"', c") is relatively 
strongly compact in W^'°° and it converge (up to a subsequence) to a solution to the original 
problem (|5.1|) . Moreover, the estimate 

sup (||/)"(t)||H={T2)+e||^^^(i)||H»(T2) + \/e||c^(t)lk=(T2)) <K 

0<t<T 

can be passed to the limit by weak lower semicontinuity and the proof is complete. D 
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Remark 5.3. The whole procedure developed in the proof of the above theorem can be easily 
generalized to the case of the third scaling introduced in section 12. 3[ 

Remark 5.4. We observe here that the power like expression for the pressure g{p) = p^ can 
be replaced by a more general one in order to achieve the same existence result as in the 
above theorem. In particular one can use g{p) = log/), thus obtaining a system which relaxes 
toward a Keller-Segel type system with linear diffusion. Therefore, some of the relaxation 
results contained in chapter H] would include Keller-Segel type system with linear diffusion 
as possible limits. This fact is not in contradiction with the finite time blow up phenomena 
occurring in the latter, because the class of initial data for which the above theorem holds is 
not significant enough in order to see the appearance of blow-up in the limit system. 
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